Strip plasmons in a two-dimensional electron gas with grounded electrodes 
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Plasmons in two-dimensional electron gas (2DEG) strips with grounded electrodes (a gate or side 
contacts) are investigated. We consider two systems: (a) the 2DEG strip with a highly conducting 
gate and (b) the 2DEG strip with semi-infinite highly conducting side contacts. The systems 
are described by hydrodynamic equations coupled with the Poisson equation. Dispersion relations 
and ac electric potential distributions are obtained. We find the plasmon modes whose potential 
distributions, and hence electron density, are localized to the edges of the strip when the wave 
number of plasmons in the direction parallel to the strip is large. Frequencies of these strip modes 
are lower than those of infinite two-dimensional plasmons. The presence of grounded electrode(s) 
significantly modifies the plasmon dispersion relations. 
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I. INTRODUCTION 

The spectrum of plasmons in electron systems depends 
on their dimensionality. In a bulk electron gas, the plas- 
mon frequency, u>, is virtually independent of the wave 
vector, q = (q x ,q y ,q z ): uj ~ yj 4Tre 2 n/mee, where e and 
m are the electron charge and effective mass, respectively, 
as is the dielectric constant, and n is the electron volume 
concentration. The spread in the electron velocities as- 
sociated with the thermal movement of electrons or the 
degeneration of the electron system results in some de- 
pendence of oj on q. However this dependence is rather 
weak pj . Stern and Chaplik Q (see also early papers 
by Ritchie Q and Ferrell discussed plasmons in the 
two-dimensional electron gas (2DEG) in the n-type in- 
version layer of metal-insulator-semiconductor structure. 
The dispersion relation for plasmons in 2DEG is given 
by @ to oc y/q. Here q = (q x ,q y ) (z-axis is directed 
perpendicular to the 2DEG), and q = \q\. Fetter 
studied plasmons in the classical 2DEG on the surface 
of liquid He. The plasmons propagating along a one- 
dimensional wire (in the 2-direction) correspond to 0, 13 
uj oc y / '\n(2/qa)q, where q = q x and a is the wire radius. 
The linear dispersion relation of 2DEG with a metallic 
gate has been found in Rcfs. HandU The first experi- 
mental observation of 2D plasmons in the 2DEG on liquid 
He was reported by Grimes and Adams Thereafter, 
Allen et al. and Tsui et al. ^2] observed far-infrared 
absorption and emission, respectively, in silicon inversion 
layers. 

Edge plasmons in 2DEG with applied magnetic field, 
so-called edge magnetplasmons, were observed by Mast 
et al. [l3L [l4| and Glattli et al. 0] Following their 
works, many theoretical investigations of edge mag- 
netplasmons in bounded structures (semi-infinite half 
plane , disk 00, and strip 00) 
were conducted using hydrodynamic electron transport 
models for 2DEGs coupled with self-consistent Poisson 
equation. The equations of these models produce integral 
equations which can be tackled by replacing their kernels 
with ones, or by expanding their solutions over complete 
sets of orthogonal functions (for comparison of these two 
methods see, for example, works by Fetter 0,0]). Edge 



plasmons in semi-infinite half plane with gatejs) were dis- 
cussed by Fetter and Nazin et al. 23] . The dis- 
cussion of early results on 2D plasmons can be found in 
the review papers by Theis 0] and Ando et al. 0] 

Rudin and Dyakonov 26] considered the edge plas- 
mons in the strip (strip plasmons) with zero magnetic 
field using a variational method, and found that two strip 
plasmon modes can exist in the strip structure (in Ref.l26l 
the lowest symmetric and antisymmetric modes were con- 
sidered). However, to use their method, it is necessary to 
choose appropriate trial functions which are usually un- 
known and hard to find a priori, especially for complex 
structures with grounded electrodes. 

In this paper, we calculate the spectra of plasmons 
propagating along strips of the 2DEG with closely lo- 
cated highly conducting grounded electrodes (a gate or 
side contacts), using the method developed in Refs. l2?l 
and l28l in which the solution of the integral equation 
reduces to the eigenvalue problem. For simplicity and 
bearing in mind 2DEG in different systems like field- 
effect transistors, we consider 2DEG without magnetic 
field assuming the uniform electron equilibrium density 
and neglecting the gradient of pressure terms in the equa- 
tion of motion [2(| • We also neglect the electron collisions 
in 2DEG, which naturally leads to some modification of 
the plasmon dispersion relation |27j and results in some 
damping. The specifics of the geometry of the 2DEG 
systems in question and the charges induced by the per- 
turbations of the electron density in the metallic region 
results in a modification of the plasmon spectra. We find 
the spatial distributions of the ac electric potential and 
the dependence of their spectra on the geometric param- 
eters of the 2DEG system. In particular, we show that 
at sufficiently large wave numbers of the plasmons, their 
ac electric potential can be localized near both edges of 
the strip. For the 2DEG strip with metallic semi-infinite 
half- plane electrodes on its both sides, we find the higher 
strip plasmon modes. 

Since the 2DEG systems similar to those studied in 
the paper can serve as resonant cavities and wave guides 
in different heterostructure devices operating in the ter- 
ahertz range of frequencies 0] , the results obtained be- 
low can be helpful of the development and optimization 
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of such devices. 



II. EQUATIONS OF THE MODEL 

We consider the following systems: (a) the 2DEG strip 
with a highly conducting gate of spacing W between 
them (gated cavity) and (b) the 2DEG strip with semi- 
infinite highly conducting side contacts (slot diode). The 
structure considered in Refs.0,|22) and l2fil can be con- 
sidered as the structure (a) with W — > oo. Schematic 
views of these structures are shown in Fig. ^ 
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FIG. 1: Schematic views of structures under consideration: 
(a) gated cavity and (b) slot diode. 

We use the hydrodynamic equations (the continuity 
equation and the Euler equation) to describe the electron 
transport: 



as 

~dt 



V • (Em) = 0, 



Ou . , e _ 
— + («• V)u = —Vcp, 
at m 



(1) 



(2) 



where £ = Y*(x,y, t) is the electron sheet concentra- 
tion in the strip, u = u(x, y, t) is the electron velocity, 
ip = ip{x,y,z,t) is the electric potential, and e and m 
are electron charge and effective mass, respectively, u is 
in xy-plane and V = (d/dx, d/dy). The hydrodynamic 
equations are coupled with the Poisson equation: 



dx 2 



dy 2 



d 2 tp iire 
dz 2 



X5(z)9(l + x)9(l-x), (3) 



where ae is the dielectric constant of the medium sur- 
rounding the strip, 6 is the Dirac delta function, and 9 is 
the Heaviside step function. 

We consider the traveling wave with frequency uj and 
wave number q y along y-axis, assuming that the potential 
(as well as other quantities) comprises both the dc and ac 
components: <p{x, y, z, t) = <fo{x, z)+ip u (x, z) exp[i(q y y— 
uit)]. Then, linearized version of Eqs. (JTJ and @ together 
with Eq. ijSJ gives the following self-consistent equation: 



d 2 ip UJ | d 2 (p, 
dx 2 
21 



+ 



dz 2 



Q v fu 



dx 2 



6{z)9{l + x)9{l-x), (4) 



where \ u — lu 2 /VI 2 and Vl — \J '2tt 2 e 2 Eo / meel is the char- 
acteristic plasma frequency. Equation Q can be reduced 
to the following: 

21 f l 

ip u (x,z) = r— / dx'Gtx, z; x', 0) 

J_( 

x (^72 -^)^(^'.O), (5) 

where G(x, z; x' , z') is the Green function which can be 
determined by the boundary conditions for each structure 
(boundary conditions and the Green function for each 
structure is introduced in subsequent sections). 

For \x\ < I and z = 0, Eq. © is reduced to the follow- 
ing integral equation: 
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1 f d 2 

dx'G(x,0;x',0) 1 

-i 



\dx' 2 



(6) 

where ^^{x) = ip LO (x,0). To solve Eq. |JB}, we expand 
ifu (x) in the cosine and sine series corresponding to sym- 
metric and antisymmetric modes of potential, respec- 
tively. The series are different for each structure because 
boundary conditions at \x\ = I are different. Then, we 
can reduce the equation to the eigenvalue problem with 
an infinite matrix whose matrix elements are determined 
by the Green function G and the expansion. The prob- 
lem can be solved numerically by truncating the matrix 
up to size TV, evaluating its elements numerically, and 
then solving the eigenvalue problem with the finite ma- 
trix. Convergence of eigenvalues and eigenvectors was 
checked by increasing N, and we use the value N = 100 
for all calculations done in the paper. An eigenvalue of 
the matrix, \ u , gives the plasma frequency of a mode 



to 
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III. STRIP PLASMONS IN A GATED CAVITY 

Since electrons in the 2DEG strip of a gated cavity 
cannot go outside, the x-component of their velocity at 
edges of the strip must be zero, i.e., in terms of the po- 
tential, dip u /dx\ z —Q\ x \—i = 0. Consequently, the Green 
function for a gated cavity is 

G ( - GC \x, z; x', z') 

= ^ K (qy^ix-x' f + iz-z 1 ) 2 ) 

-j-K (q y ^{x - x') 2 + (z-z'- 2W) 2 ) , (8) 

where the first term of Eq. © is the Green function for 
a free cavity. The expansions which satisfy the boundary 
conditions are the following: 



¥>L s) (z) 



V2 



k=l 



k cos(q 2k x) 



(9) 
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for symmetric modes and 

oo 

^\x) = Y J c { £ ) sHd2 k -ix) (10) 

fc=l 

for antisymmetric modes, where q n = nn/2l. The ef- 
fective wave number corresponding to the x-direction is 
equal to q2n for n-th symmetric mode (as a special case, 
zero-th mode correspond to q 2n — 0) and to tfen-i for n- 
th antisymmetric mode. Substituting Eqs. © and 11U|) 
into Eq. JfjJ) and using the orthonormality of each term 
in expansions, we arrive at the following eigenvalue prob- 
lem: 

oo 

E W = WW en) 

fe'=0 

for symmetric modes, where fc = 0, 1, 2, • • • and 

f# = \{lyl? j\d£ J^dCG^^a (12) 
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FIG. 2: Dispersion relation for a free cavity (W/l — > oo) 
with <j„//7r = 0, 1/2, 1, 3/2, 2, 5/2. Solid lines with lowest 
frequencies through highest ones correspond to q n l/n — 
through q n l/Tr — 5/2. The dotted line is the dispersion rela- 
tion of infinite 2D plasmons. Two modes whose frequencies 
are lower than those of infinite 2D plasmons for g s l > 1 are 
identified as strip plasmon modes. 



x cos(q 2k lO (13) 

when fc > 0, 

tfffi - — [M 2 + (%o 2 ] ^ r tfo^&o 

x cosfaylZ') (14) 
when /c' > 0, and 

x cos(q 2k l0 cos(q 2 k>l£) (15) 
when fc, fc' > 0. For antisymmetric modes, we obtain 

oo 

E^^cf, (16) 
fc'=l 

where fc = 1, 2, 3, • • • and 

~ J-i j-i 

x sin(g 2 fc-i/0 sin(g 2 fc'-i^')- (I 7 ) 

Here we use the notation G( GC )(£, £') = 
G (GC )(Z£,0;Z£',0). The details of numerical com- 
putation of Eqs. (|12|) - (|15fl and H17|) are discussed in the 
Appendix. 

Figure |21 shows the dispersion relation calculated for a 
free cavity. The dotted line is the dispersion relation of 
infinite 2D plasmons which is given by lo p = W q y l/i:fl. 



Frequencies in the limit q y l — > coincide with those cal- 
culated for standing plasma waves in a free cavity in 
Ref. |2^. It can be seen from Fig. [21 that the lowest sym- 
metric and antisymmetric modes have lower frequencies 
than those of infinite 2D plasmons for q y l^> 1, while fre- 
quencies of higher modes approach to those of infinite 2D 
plasmons. These two modes are identified as strip plas- 
mons, as is also evident from their potential distributions 
(Fig. I3- This phenomenon was already discussed theo- 
retically in the case of a semi-infinite half-plane |18| and 
in the case of a free cavity (2||. For q y l~^> 1, the effect of 
one edge on another can be neglected, and the difference 
between symmetric and antisymmetric modes becomes 
unimportant. The lowest modes in our case are equiva- 
lent to an edge mode in the case of a half-plane. The dis- 
persion curves for lowest symmetric and antisymmetric 
modes are identical with those obtained in Ref. |2£j The 
higher modes seem to correspond to the normal wave-like 
modes in x-direction whose frequencies are proportional 

t° \Jln + 1y because the effective wave number q n can 
be considered to take continuous value when q y l 1. 

The dispersion relation for plasmons in a gated cav- 
ity is illustrated in Fig. 21 For the zero-th symmetric 
mode (q n l/n = 0), one can see that, as W/l decreases, 
the plasma frequencies decrease and their dependence on 
q y l becomes linear. When q y l/W, the effect of the 
gate is small and the frequencies approache to those for 
a free cavity (see curves for W/l = 0.1). For the first 
symmetric mode (q n l/^ = 1) the dependence is more 
complicated due to nonzero frequency of plasmons in the 
limit q y l — > and its dependence of W/l [2q . Figure El 
shows the plasma frequencies as a function of W/l for 
fixed q y l. For q y l = 1, the influence of the gate on the fre- 
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FIG. 3: Potential distributions of lowest antisymmetric mode 
for a free cavity with q y l = 1 (above) and q y l = 10 (below). 
The localization of potential at the edges of the 2DEG strip 
becomes apparent for q y l S> 1. 




W/l 

FIG. 5: Plasma frequency vs the ratio W/l with q y l — 1, 10 
and different q„l/n for a gated cavity. For each q y l, four solid 
lines with lowest frequenciy through highest ones correspond 
to q n l/^ = 0, 1/2, 1, 3/2. For q y l = 10, curves for q„l/n — 
and q n l/n = 1/2 are indistinguishable. 



IV. EDGE PLASMONS IN A SLOT DIODE 



quencies is gradual, i.e., the frequencies slowly decrease 
with decreasing W/l and tend to zero as W/l — > 0. On 
the contrary, for q y l = 10, the change in the frequency 
from zero to the value corresponding to W/l — > oo is 
more rapid. It is also worth mentioning that for q y l = 1 
the frequency differences between lowest and second low- 
est symmetric and antisymmetric modes are significant 
whereas they are not the case for q y l — 10. 




FIG. 4: Dispersion relation for a gated cavity with q n l/n — 
0, 1 and W/l = 0.01, 0.025, 0.05, 0.1, oo. 



For a slot diode, the boundary conditions are given by 
fuj\\x\>i,z=o = 0; an d the relevant Green function is 

G (SD) (x,z;x',z') 
= ^ E K o (q y ly/(O-0'-2mr) 2 + ty-il/) 2 ) 

71 — — OO 
1 OO 

E K ° (y V(0 + W- 2wrp + - #) 2 ) , 

Tl — — OO 

(18) 

where 0+iip — cos^^ir+iz)//] and 9'+iip' — cos~ 1 [(a;' + 
iz') /I] . When evaluating Eq. (|18fl numerically, series are 
cut at appropriately large n, so that the ra-th term can 
be neglected. Using the expansion 



<PuH x ) = E c fe cos fe,2fc-ia; 



k=l 

for symmetric modes and 

oo 



<PwH x ) = E c fe sin fe,2fca;) 
fc=i 



(19) 



(20) 



for antisymmetric modes, we arrive at the eigenvalue 
problem similar to Eqs. i|ll|) and (|16|) [only a difference 
is that the indices k and k' start from 1 in Eq. Ijllfl ] with 
the matrix elements 

tfffi = -[(qw-ilf + iqyl) 2 ] f f dt'GW&O 
~ J -i J -i 

x cos(q 2 k-ilO cos (<l2k<-il£') (21) 
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for symmetric modes and 



g(a) _ 
7 kk' — 



-[(q x ,2k'l) 2 + (q y l) 2 ] / d£ 
n J -l 

x sin(<fe.,2fcZf) sin(g X)2 fe'^')- 



(22) 



Solving the eigenvalue problems numerically and sub- 
stituting eigenvalues into Eq. Q, dispersion relation for 
a slot diode is obtained (Fig. |HJ). Frequencies in the 
limit q y l — > coincide with those calculated for stand- 
ing plasma waves in a slot diode in Ref. |2^. Aside from 
the first modes, we have found that higher strip modes 
can exist. By observing the potential distributions, it 
was turned out that the potential distribution of a mode 
begins to change and be localized near the edges of the 
2DEG strip at the wave number q y at which the disper- 
sion relation curve of the mode crosses over the curve 
of infinite 2D plasmons. One can also see from Fig. 
that dispersion relation curves of symmetric and anti- 
symmetric modes with same index begin to coincide at 
this wave number. Figure0shows potential distributions 
of symmetric modes for q n l/"K = 1/2,3/2, and 5/2 with 
q y l = 30. From Fig. [7] one can see the increasing of peak 
numbers as qnl/ir increases, and peaks become larger and 
wider as they are farther from an edge of the 2DEG strip. 





FIG. 7: Potential distributions of three symmetric modes 
{q n l/n = 1/2,3/2, and 5/2.) with q y l = 30. These strip 
modes have two, four, and six peaks, respectively. 



the spatial distributions of ac electric potential. In par- 
ticular, it has been shown that at sufficiently large wave 
numbers of plasmons, their ac electric potential of some 
modes can be localized near both edges of the strip, and 
their frequencies are lower than those of infinite 2D plas- 
mons. 
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q I 

FIG. 6: Dispersion relation for a slot diode with q n l/i: = 
1/2, 1, 3/2, 2, 5/2, 3. Solid lines with lowest frequencies 
through highest ones correspond to q n l/n = 1/2 through 
q„l/n — 3. The dotted line is the dispersion relation of infinite 
2D plasmons. 



V. CONCLUSIONS 

In summary, we have studied plasmons in the following 
2D systems: (a) the 2DEG strip with a highly conducting 
gate (gated cavity) and (b) the 2DEG strip with semi- 
infinite highly conducting side contacts (slot diode). We 
have found the dispersion relations of the plasmons and 



Eqs. H12[) through (|14f> can be reduced to the following 
simpler expressions which contain only single integrals: 

*oo = -ily 1 ) 2 {K (2q y l)C- 1 {2q y l) + ^ 1 (2g„Z)£ (2q„/) 
- [1 - 2q y lK 1 {2q y l)} /2Tr(q y l) 2 

-^/ 2 ^( 1 -|)^o(^v / ?+^)} ! (AT) 

where K n is the modified Bessel function of second kind, 
C n is the modified Struve function and a = 2W/1, 

= (-l) k+1 ^(q y l) 2 £ G( GC \0sHk*0, (A.2) 

and 

= (-D fe,+1 ^[(^) 2 + (%o 2 ] 

r 2 

x / G (GC) (£)sin(fcV0. (A.3) 
Jo 
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Here, G<- Gc ^{£) = G (GC )(£,0). Eq. JTSJ can be reduced 
to 



and 



w = (-i)^^ 4[(fc ;f^f ] [T fc (2)-r,(2)] 

(A.5) 

^ (a - 6) 



bV) 2 + ( gj/ n2 

(fc 2 - k' 2 )n 

when fc 7^ fc'. Similarly, Eq. 117|) can be reduced to 



and 

g cg = (-d^i W + faf ] [Tft(2) _ Tfe , (2)] , 

(A.7) 

where ft = (2k- l)/2 and ft' = (2fc' - l)/2, when fc 7^ fc'. 
Here, 

5,(0 = jT 4 # sin(z,<') A G (GC) (e0 (A< 



and 



T„(0 = £ de'Icos^O - 1]^G (GC) (£')- (A.9) 



Since all integrands of integrations in Eqs. I|A.1|) through 
(|A.9|I have no singularity, the Gaussian quadrature j3j| 
was used to perform the numerical integrations. 

On the other hand, since Eq. (|18|1 is very complicated, 
we directly performed double integrations in Eqs. 121|) 
and H22f> after integrating them by parts: 



it = [fa-xl? + (lylf] 

X I d ^' g(SD) ^>OMQ2k-il0^s(q 2k/ - 1 ie) 

(A.10) 



and 



1 <% [ deg (SD \z,Ocos(q 2k i0sm(q2k'ie), 

-1 J-l 

(A.ll) 



where G<- SD ^ {£,,£) = G<- SD \t,C)dt. Since G< SZ3 > has 

logarithmic singularity at £ = £', Q^ SD "> is the continuous 
functioned £ and £'. g( SD ) was calculated using adaptive 
quadrature |3l|]. Finally, Eqs. (|A.10() and l|A.lll) were 
calculated by Simpson's double integral method [3l| . 
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